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MARTIN KERNELS FOR MARKOV PROCESSES WITH JUMPS 


TOMASZ JUSZCZYSZYN, MATEUSZ KWASNICKI 


Abstract. We prove existence of boundary limits of ratios of positive harmonic func¬ 
tions for a wide class of Markov processes with jumps and irregular domains, in the 
context of general metric measure spaces. As a corollary, we prove uniqueness of the 
Martin kernel at each boundary point, that is, we identify the Martin boundary with 
the topological boundary. We also prove a Martin representation theorem for harmonic 
functions. Examples covered by our results include: strictly stable Levy processes in R d 
with positive continuous density of the Levy measure; stable-like processes in R rf and 
in domains; and stable-like subordinate diffusions in metric measure spaces. 


1. Introduction 


The purpose of this article is to study boundary limits of ratios of positive functions 
which are harmonic in an arbitrary open set with respect to a Markov process with jumps. 
The proof of our main result, Theorem 13.11 relies on the boundary Harnack inequality 
for Markov processes with jumps, proved recently in |12| . and the oscillation reduction 
argument, developed in [5] and |llj . As an application, we obtain Martin representation 
of harmonic functions in Theorem 13.61 

To explain the motivation for our research, we begin with a discussion of the classical 
case, where harmonicity has its usual meaning: / is harmonic in an open set D if A/ = 0 
in D. The boundary Harnack inequality is a statement about positive harmonic functions 
in an open set, which are equal to zero on a part of the boundary. The result states that 
if D is regular enough (for example, a Lipschitz domain), Xq is a boundary point of D, f 
and g are positive and harmonic in D , and both / and g converge to 0 on dD C\B{x o, R), 
then for every r G (0, R) the ratio f/g has bounded relative oscillation in D D B(xo,r): 
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Here c = c(D,xo,r, R) is a constant that depends only on the local geometric properties 
of D near x 0 , and B(xo,r) denotes the ball of radius r, centred at Xq. The boundary 
Harnack inequality was first proved independently by A. Ancona ([5]), B. Dahlberg ([IB]) 
and J.-M. Wu ([3T|) for Lipschitz domains, and then extended by numerous authors to 
a wider class of domains and elliptic operators. We refer to in El El SI m for further 
discussion and references. 

Under appropriate assumptions on the regularity of D , the estimate (11.11) turns out to 
be self-improving as r —» 0 + , in the sense that the constant c in (11.11) converges to 1 as 
r — y 0 + . Equivalently, the boundary limit 


lim 

x—>xo 
xeD 


f(x) 

9{x) 


( 1 . 2 ) 


Date: September 21, 2015. 

Key words and phrases. Markov process, jump process, killed process, boundary Harnack inequality, 
boundary limit, Martin kernel, Martin boundary, Martin representation. 

Work supported by the Polish National Science Centre (NCN) grant no. 2011/03/D/ST1/00311. 

1 









2 


TOMASZ JUSZCZYSZYN, MATEUSZ KWASNICKI 


exists. When D is a Lipschitz domain, then in fact c(D , xq, r, R ) is of order as r —* 0 + 
for some j3 > 0, which means that f/g extends to a Holder continuous function at Xq. 

A closely related concept of Martin representation of positive harmonic functions was 
introduced by R. S. Martin in his beautiful article [28], more than two decades before the 
boundary Harnack inequality became available. Given the existence of limits (11.21) (for 
example, if D is a Lipschitz domain), Martin’s result asserts that there is a one-to-one 
correspondence between positive harmonic functions / in D and positive measures /i on 
the boundary of D. The two objects are linked by the formula 


f(x)= / M D (x,z)n(dz), 

JdD 

where the Martin kernel is defined as the boundary limit of the ratio of Green functions : 

G D (x,y) 


Adr>(x,z) = lim 


K c " [i - 


(1.3) 


Here x G D is an arbitrarily fixed reference point. 

One of numerous equivalent definitions of harmonicity links harmonic functions with 
the Brownian motion: / is harmonic in D if and only if / has the mean-value property 
with respect to the distributions of the Brownian motion X t at first exit times: 


f(x) = E x f(X(ru)) (1.4) 

for all bounded open sets U such that the closure of U is contained in D. Here E., ; denotes 
the expectation (and P x will denote the probability) corresponding to the Brownian 
motion process X t that starts at x, and Tu is the time of first exit from U: 


T V = inf {t > 0 :X t <£U}. 

This probabilistic definition has a number of advantages: it extends immediately to 
general Markov processes X t , and it captures easily boundary conditions imposed on 
harmonic functions. More precisely, in the general statement of the boundary Harnack 
inequality one requires that positive harmonic functions / and g converge to zero at 
each boundary point in dD fl B(;Xq^R) that is regular for the Dirichlet. problem. This 
condition translates to requiring that (11.41) holds for all bounded open sets U such that 
U C D U (dD D B(xq,R)), with no reference to the notion of regular boundary points. 
Here we understand that f — g — 0 in dD fl B(x 0 , R). 

In this article we are interested in Markov processes with jumps, and from now on 
by saying that a function is harmonic we understand that it has the mean-value prop¬ 
erty ( 11 . 41 ) with respect to a Markov process X t with jumps. In this case in order to 
evaluate f(X(ru )) in ( 11 . 41 ) the function / needs to be defined everywhere, not just in D. 
For this reason one needs to replace the boundary condition / = g — 0 in dDnB(xo, R) in 
the statement of the boundary Harnack inequality with the exterior condition f = q = 0 
in D c fl B(xq, R). 

The history of the boundary Harnack inequality for Markov processes with jumps 
starts with the article by K. Bogdan ([6]), where he proved the result for the isotropic 
stable Levy process (equivalently: for the fractional Laplace operator — (—A) Q//2 ) anc l 
Lipschitz domains. Later this was extended to more general sets ( |301 IT] ) and processes 
([HI m m |2ll [22j [23l [2U [25]). Recently, a rather general result for Markov processes 
with jumps was proved in [T2J, and this is our starting point in the study of boundary 
limits ( 11 . 21 ) . 

The existence of the boundary limit ( 11 . 21 ) in this context was first proved independently 
by K. Bogdan (0) and by Z.-Q. Chen and R. Song (H) for the isotropic stable Levy 
process and Lipschitz domains. This required an appropriate modification of the classical 
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reasoning due to the presence of jumps. Since then essentially every time the boundary 
Harnack inequality was established for a given Markov process with jumps in a given 
class of domains, the existence of boundary limits (11.21) followed; see |26| for the most 
recent result of this kind. With one exception, however, the class of open sets under 
consideration was always limited to certain disconnected analogues of non-tangentially 
accessible domains, typically called fat sets. The single more general result is proved 
in m for the isotropic stable Levy process, where completely arbitrary open sets are 
allowed. 

For the existence of boundary limits, we follow the approach of HH using the boundary 
Harnack inequality of [T2], and prove in our main results, Theorems 13.11 and 13.61 the 
existence of boundary limits of ratios of harmonic functions for arbitrary open sets and 
rather general Markov processes with jumps, as well as Martin representation of such 
functions. The application of the method developed in m in the present setting requires 
significant modifications. Further changes are introduced in order to make the description 
of the proof more accessible; for example, we first give a simpler argument which does 
not assert uniform convergence with respect to the domain of harmonicity, and only then 
explain how one improves it to get a domain-uniform version. 

The proof of Martin representation theorem for the isotropic stable Levy processes 
in m is self-contained. It is possible to extend the method of m to our general setting, 
but that would require rather lengthy and technical arguments. For this reason, unlike 
in HU, we refer to the general theory of Martin boundary. Our argument still requires 
extension of some elements of HU for more general Markov processes, but the most 
involved part of the proof is avoided. For an excellent exposition of the general theory of 
Martin boundary, we refer to Chapter 14 of pg. 

We conclude the introduction with a description of the structure of this article. The 
assumptions for the boundary Harnack inequality of j 1 2] are briefly recalled in Section [21 
We omit a detailed discussion of these conditions and refer the interested reader to the 
original paper. Instead, we present a number of examples right after the statement of 
Theorems 13.11 and 13 .6 1 in Section [3] We also provide a counter-example, which shows that 
the boundary limits (11.21) typically fail to exist in irregular domains when the process X t 
has a noil-trivial diffusion part. Finally, in Section [4] we prove Theorems 13.11 and 13.61 


2. Fundamental assumptions for the boundary Harnack inequality 

The formal statement of the assumptions for Theorem 13.11 requires some effort. We 
assume that (X, d, m) is a locally compact metric measure space in which all bounded 
closed sets are compact and m has full support, and that Rq > 0 (possibly R 0 = oo) is a 
localisation radius such that X \ B(x, r) ^ 0 if x e X and 0 < r < 2 R 0 . 

In p2] the following four conditions are introduced. A detailed discussion of these 
assumptions is beyond the scope of the present article, we refer the reader to pj2] for 
more information. Here we only state the conditions, without explaining in a formal 
way the notions of semi-polar and polar sets, processes in duality X t and X tl their 
generators 21 and 21, densities v(x, y ) and u(x, y) (with respect to the measure m ) of the 
Levy kernels of X t and X t , as well as their Green functions Go^x^y) = GD(y,x). We 
note that v(x,y) describes the intensity of jumps from x to y and it is commonly used 
throughout the article. The Green function Gn^x^y) is required for Theorem 13.61 only: 
informally, Gd(x , y) is the average amount of time spent near y by the process X t , started 
at x , until t d . 
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Assumption A. Hunt processes X t and X t are dual with respect to the measure m. The 
transition semigroups of X t and X t are both Feller and strong Feller. Every semi-polar 
set of X t is polar. 

Assumption B. There is a linear subspace D of 2D (21) D 2) (21) satisfying the following 
condition. If K is compact, D is open, and K C D C X, then there is / G 2D such that 
f(x) = 1 for x G K, f(x) = 0 for x G X \ D, 0 < f(x) < 1 for x G X, and the boundary 
of the set {a: : f(x)> 0} has measure m zero. 

Assumption C. We have v{x,y) = 0(y,x) > 0 for all x, y G X, x ^ y. If x 0 G X , 
0 < r < R < R 0 , x G B{x 0 , r) and y G X \ B(x 0 , R ), then 

C Uvy U ( X 0> y) ^ U ( X > y) ^ C LevyV(x 0 , y), C^v{x 0 , y) < f>(x, y) < C' L evyh(.X'o, y), (2.1) 

with C^Levy f-^Levy (^Oj ^ . 

Assumption D. If xo G X, 0 < r < s < R < i? 0 and B = B(x o, R), then 

C gTeen = C green (x 0 ,r,s,R) = sup sup max(G B (x, y), G B (x, y)) < oo. (2.2) 

x£B(xo,r) y£X\B(xQ,s) 

We denote 

p(K, D) = inf sup max(2l/(x), 2l/(x)), (2.3) 

/ x&X 

where the inhmum is taken over all functions / described by the Assumption iBl If xq G X 
and 0 < r < R < R 0 , then we denote 

CLevy-infOro, r, R) = _ inf mm(u(x 0 , y ), u(x 0 , y )), 

yeB(xo,R)\B(xo,r) 

and 

C exit (x 0 ,r) = sup max(E x r B ( xo ,r)j ®a:TB(xOjr)). 

x£B(xo,r) 

Following j6], we say that / is a regular harmonic function in an open set D if the mean- 
value property (11.41) holds with U = D. By the strong Markov property, this implies 
that (11.41) holds for arbitrary open U C D, so in particular / is harmonic in D. The 
following theorem is a reformulation of the main result of |12j . 

Theorem 2.1 (Lemma 3.2 and Theorems 3.4 and 3.5 in [T2]). Suppose that xo G X, 
0 < ri < r 2 < r 3 < r$ < Rq and a non-negative function f is a regular harmonic function 
in D fl B(xo,r 6 ), which is equal to zero in B(xo,re) \ D. Then 


f{x) ~ CBHlE x T£>nB(:E0,r2) 


f(y)u(x 0 ,y)m(dy) 


'X\B(x 0 ,r 3 ) 

for x G D fl B(xq, rf), where Cbhi = Cbhi^o, f l, r 2, r 6 ) is defined as 

Cbhi = C Uvy (xo,r 2 ,r 3 ) + 2 p(B(x 0 ,r 3 ) \ B(x 0 , r 2 ) , B (x 0 , r 8 ) \ B(x 0 ,r 1 ))x 


x G green (x 0 ,r 3 ,r 4 ,r 6 ) + 


Cexit (^0 j ) ( C*Levy (*^0 j ^4 j ^5 ) ) 

m(B(x 0 ,u)) 

f p(B(x 0 ,r 5 ),B(x 0 ,r 6 )) 

x —y-7-7- h C h ^ V y{xo,re,r 7 )m(B(x 0 ,r 6 )) 

V F L6v y-inf(.x'o, r 5 , r 7 ) 

for some r 4 , r$, r 7 , r 8 such that 0 < r\ < r 2 < r 3 < r4 < r 5 < r 6 < r 7 < r 8 . 


x 
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Note that it is important that / is non-negative everywhere, not just in D. Theorem 12.II 
implies the more classical statement of the boundary Harnack inequality (Theorem 3.5 
in pj2J): if / and g satisfy the assumptions of Theorem 12.11 then 


sup 

X^D r 


f{x) 

g{%) 


— Cbhi i n f 

xeD 7 


f(x) 

a(x) 


(2.4) 


as in 01. ip . We remark that although the original statement allows for an arbitrary 
sequence of radii, it will be sufficient for us to consider rr = r, r 2 = 2r, r% = 3r and 
r 6 = 4r, and we will commonly write Cbhi = Cbhi(to,f) = Cbhi(xoj t, 2r, 3r, 4r) in this 
case. 


3. Main results and examples 


For the existence of limits, we introduce one more definition. If Xq G X and 0 < r < 
R < Rq, we let 


Cbevy-int C^Levy-int (^ 0 1 T *5 -R) 


Jx\B{ X0 ,r) u ( x ^y) m ( d y) 
Ix\B( X0 , R ) u ( x ^y) m ( d y)' 


(3.1) 


We use a short-hand notation f ~ eg for the two inequalities c l g < f < eg, where c > 0 
is a positive constant. 


Theorem 3.1. Let D C £ be open, xq G dD and R > 0. Suppose that: 

(i) X t satisfies Assumptions HI through [D|: 

(ii) lim C L6vy (x 0 ,r,i?) = 1; 

r—>■ 0+ 

(iii) the constant Cjj svy (xo,r, 2 r) is bounded in r, 0 < 2 r < R 0 ; 

(iv) the constant C^ vy .- mt (x 0 ,r,2r) is bounded in r, 0 < 2r < R 0 ; 

(v) the constant Cbhi(^o, r , 2r, 3r, 4r) is bounded in r, 0 < 4r < R 0 . 

Suppose furthermore that non-negative functions f and g are regular harmonic functions 
in DnBfx o, R) and are equal to zero in B(x o, R)\D. Then either one of f and g is zero 
everywhere in D, or the finite, positive boundary limit of f(x)/g(x) exists as x Xq, 
x G D. Furthermore, 


Hm /(x) = J X \B(o,r) u ( x ^y)f(y) m ( d y) 
XGD 3(' X ) ^0+ fx\B( 0,r) U ( X 0> y)9(y) m (dy ) ’ 


(3.2) 


Remark 3.2. Condition 
acterised by the property 


is required only for inaccessible boundary points Xq, char- 
E y TDnB(x 0 ,R) m ( d y) < oo. The result for accessible 

and 


' DnB(xo,R) 

boundary points Xq, for which the integral is infinite, holds under conditions 
through 


m 


Remark 3.3. Theorem 13.11 also holds with g(x) = Ri x TDnB(x 0 ,R)- This is formally shown 
in Section 14.41 but the informal explanation is rather straightforward: g is essentially a 
regular harmonic function in D C\B(x o, R) (in sharp contrast with the case of continuous 
Markov processes). 

Indeed, suppose that X is unbounded, D is a bounded open set and that Chevy^o, r , R) 
converges to 1 as R —> oo. By Dynkin’s formula (see Lemma l4~2l and estimate (14.21) below), 

F . = ,. P.(A :(T D )£X\B(z a ,R)) 

D U\B(x 

is the limit of regular harmonic functions in D. Since the estimates in Theorem 13.11 are 
uniform in / and g, we obtain the desired result. (Note that the formal argument is 
completely different and requires no further assumptions on X and X t .) 
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Remark 3.4. As remarked in the introduction, the limit in (13.21) exists if and only if the 
relative oscillation of / and g converges to one, that is, 

Um SUPxgDnB(x 0 ,r)(/(^)M^)) _ 
r^o+ inf x€DnB (x 0 , r )(f(x)/g(x)) 

By inspecting the proof of Theorem 13 .1 1 one immediately sees that, given D and Xq, the 
boundary limits exist uniformly in / and g, in the sense that 

sup a , e r,nB(xo,r)(/(®)/p(®)) 


lim sup 

r ^ 0+ f,g mi xeDnB(x 0 ,r){f{x)/g{x)) 


= 1, 


with the supremum taken over all / and g satisfying the assumptions of the theorem. We 
remark that in fact one can prove uniformity also in D , just as in mi, by appropriately 
modifying the final part of the proof. More formally, 


sup xeD nB(x 0 ,r)(f(x)/g(x)) , 

lim sup ——- , - , , —— = 1, 


(3.3) 


>-W + DJ.g in!x€DnB(i 0 ,r)(f(x)/g{x)) 

where the supremum is taken over all open sets D and / and g satisfying the assumptions 
of the theorem (here we let the ratio sup / inf be equal to 1 if D fl B(xo,r) is empty). 
The proof of this result is sketched in Section 14.41 


Remark 3.5. It is not necessary to assume that Xq G 8D in Theorem 13.11 For x () f I) 
the statement is void, but for Xq G D we obtain relative continuity of positive harmonic 
functions: if / and g are positive harmonic functions in D , then f/g is continuous in D. 
If the process is conservative, then the constant 1 is harmonic, and consequently positive 
harmonic functions are continuous. If in addition the characteristics of the process (that 
the constants in conditions (ii) through (v)) do not depend on xq, then harmonic 


is, 


functions are in fact uniformly continuous (see also Remark 13.41) . 


Before we discuss examples, we provide one application. Recall that the Green function 
Gn(x,y ) is the density of the mean occupation measure of X t up to td, that is, 


rr D 


G D (x,y)m(dy) = E x l A (X s )ds. 


J A Jo 

Under Assumptions lAl and iDl there is a version of Gf)(x,y) which is a harmonic function 
of x G D \ {y}, and a co-harmonic (that is, harmonic for the dual process) function of 
y G D \ {x}. Hence, Theorem 13.11 (or, more precisely, its version for the dual process) 
immediately implies the existence of the Martin kernel 


M d (x,z) 


,. G d (x, y) 

lim —- . 

G D (x,y) 


for z = x o (this is exactly the same as the classical definition (ll.3jl ). Informally, the Mar¬ 
tin boundary 8mD of a set D is the set of all possible ways a point y G D approaches the 
boundary in such a way that the ratio Gd(x, y)/Gc{x, y) converges for every iGh (with 
arbitrarily fixed x G D). More formally, D U DmD is the Constantinescu-Cornea com- 
pactification of D with respect to the family of functions {Gd(x, -)/Gd(x, •) : x G D}: the 
smallest compact space which contains D and on which these functions have continuous 
extensions. 


Theorem 3.6. Let D C X be bounded and open, and if X is compact, then assume in 
addition that E x (tjj) is finite and bounded in x G D. Suppose that the assumptions of 
Theorem \3.1\ are satisfied uniformly for all Xq G D. Then the following assertions hold. 
(a) The Martin boundary 8 mD coincides with the topological boundary 8D. 
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(b) The Martin kernel M B fx, z ) is a harmonic function in D with respect to x if and 
only if z is an accessible boundary point: f DnB ^ XQ U) ^ y ' r DoB{x 0 ,R)i ri '{dy) = oo. 

(c) In this case M B {x, z) is a minimal harmonic function: if f is a harmonic function 
in D and 0 < f(x) < Md{x, z ) for all x G X, then f(x) is a multiple of M B (x, z ). 

(d) Every non-negative function f which is a harmonic function in D has a unique 
representation 

f{x)= I ( [ G D {x,y)u(y,z)m(dy ) 

JX\(DUdmD) \Jd 

+ / M D (x, z)n(dz), 

JdmD 

where p is a measure on d m D, the set of accessible boundary points of D. 

(e) Conversely, given any non-negative function f and any measure p on d m D, the 
right-hand side of (13.411 is either a harmonic function in D or infinity everywhere 
in D. 


f(z)m(dz) 


(3.4) 


Remark 3.7. The terms accessible and inaccessible correspond to the probabilistic the¬ 
ory of Martin boundary. To be specific, the process X t killed at the time of first exit 
from D and conditioned in the sense of Doob by the Martin kernel M B (-, z) converges at 
its lifetime to 2 when z is accessible, and dies out in D when z is inaccessible. We refer 
to ra for more information. 

Remark 3.8. Unlike in the case of isotropic stable Levy processes in m, description of 
the infinite part of the Martin boundary of D for unbounded open sets is a completely 
different problem. 


The boundary Harnack inequality stated in Theorem 12.11 was applied to a variety of 
Markov processes in Section 5 of [12]. The scale-invariant version of Theorem 12.11 under 
a-stable-like scaling discussed therein already asserts conditions 
rein 13.11 Verification of the remaining conditions 


n 


m 


and (v) in Theo- 


and (iv) is typically straightforward, 


and we obtain several classes of processes for which Theorems 13.11 and 13.61 apply. 

In our the first example, we use the result of Example 5.5 in [12] . where boundary 
Harnack inequality for Levy processes is considered. In the asymmetric case, equality 
of the notions of semi-polar and polar sets (in Assumption 0 is not trivial, and this 
was apparently overlooked in m Fortunately, for all asymmetric Levy processes listed 
therein, this condition is satisfied by Theorem 2 in [29] , 


Example 3.9 (Strictly stable Levy processes). Let m be the Lebesgue measure in R^, 
Ro = oo. Suppose that X t is a strictly a-stable Levy process in R d , where d > 1 and 
0 < ol < 2. Suppose, furthermore, that the Levy measure of X t has a density function of 
the form u(z) = (p{z/\z\)\z\~ d ~ a , with <p continuous and positive on the unit sphere (for 
Levy processes, v(x,y) = v(y — x)). It is easy to see that C^ vy (x 0 , r, R) converges to 1 
as r —» 0 + and that C^ vy . int (x 0 ,r, R) = ( R/r) a . By Example 5.5 in [T2], X t satisfies the 
other assumptions of Theorem 13.11 and so we may use Theorems 13.11 and 13.61 


We remark that the above example can be extended to more general Levy processes, 
including many subordinate Brownian motions and, more generally, unimodal isotropic 
Levy processes. This is based on estimates obtained recently in |Q| HQ, [T7J HQ] and will be 
studied in detail in |T8]. Another extensions can be obtained by allowing the Levy kernel 
to depend on x or restricting it to a domain, as described in the following two examples. 
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Example 3.10 (Stable-like processes). Let m be the Lebesgue measure in B d , R 0 = oo. 
Suppose that 0 < a < 2 and 

u(x,y) = p(x,y)\x - y\~ d ~ a , 

where <p is symmetric (that is, <p(x, y) = ip(y, x)), bounded by positive constants, smooth, 
and has bounded partial derivatives of all orders. As in Example 5.6 in ns, in this case 
there is a pure-jump process X t with the Levy kernel is(x,y)m(dy), and the assumptions 
of Theorem 13.11 are satisfied. 


Example 3.11 (Reflected stable processes). Let 0 < a < 2. Let X be the closure of 
either a Lipschitz domain in R fl if a < 1 or a C l,a+£ domain in R (/ if a > 1 (with some 
e > 0). Let m be the Lebesgue measure on X, and v(x, y ) = c\x — y\^ d ~ a for some c > 0. 
Again as in Example 5.6 in da. there is a pure-jump process X t with the Levy kernel 
v(x,y)m(dy), and the assumptions of Theorem 13.11 are satisfied for some Rq > 0. 

The state space X need not be Euclidean. 


Example 3.12 (Stable-like subordinate diffusions). Let X be a sufficiently regular metric 
measure space in which there exists a diffusion process. For a rigorous definition, we refer 
to Example 5.7 in m examples include Riemannian manifolds, Sierpihski gaskets or the 
Sierpihski carpet. Suppose that 0 < a < d w , where d w is the walk dimension of X (that 
is, an approximate scaling exponent for the diffusion process). Finally, let X t be a process 
subordinate to the diffusion process, corresponding to the (a/d w )-stable subordinator. In 
Example 5.7 in [T2] it is shown that X t satisfies conditions (i), (iii) and (v) of Theorem 13.11 
and one easily proves that Chevy-int^o, A R) < c(R/r) a for some c > 0. Verification of (ii) 
requires some work, especially when X is unbounded. For this reason, we only sketch the 
argument for compact X. For some c > 0 we have 

poo 

v(x,y) = c t~ 1 - a/dv q t (x,y)dt, 


where q t (x,y) is the transition density of the diffusion process. Since for each t > 0, q t is 
Holder continuous, it is easy to see that v(x, y) is positive and uniformly continuous in 
x G B(x 0 ,r), y G X \ B(xq,R), which clearly implies condition (ii). It follows that The¬ 


orems 13.11 and 13.61 apply to stable-like subordinate diffusions in compact metric measure 
spaces. 


Surprisingly, Theorem 13.11 is not influenced by killing. 


Example 3.13 (Processes with a multiplicative functional). Let M t be a strong continu¬ 
ous multiplicative functional such that M 0 = 1 with probability one for all starting points 
x 6 X. Such a functional describes gradual killing of the process X t , and is typically ob¬ 
tained as the Feynman-Kac functional M t = exp(— V(X s )ds) for some non-negative 
function V. A function / is said to be harmonic with respect to the pair (Xt, M t ) if it 
has the mean-value property 

f(x) = ~E x (f (X (t v ))M ( Tu )) 

instead of (II.4[) . As in Theorem 5.10 in [12] . if the assumptions of Theorem 13.11 are 
satisfied by the process X t , then the conclusion also holds for functions harmonic with 
respect to the pair (X t ,M t ). 


Our final example shows that when X t has non-vanishing diffusion part, one cannot 
expect the existence of boundary limits (11.21) unless some geometric restrictions on D are 
imposed. For corresponding positive results in smooth domains, see [23] . 
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Example 3.14 (Mixture of Brownian motion and stable process). Let X = R and let m 
be the Lebesgue measure. Let X t be a one-dimensional Levy process which is the sum 
of two independent Levy processes: the Brownian motion and the symmetric a-stable 
Levy process for some a G (1,2). That is, the characteristic exponent of X t is given by 
Ci £ 2 + C 2 |£|" for some ci, C 2 > 0. Denote D = (—1,1)\{0}. Let p t (y — x ) be the continuous 
version of the transition density of X t . Then the three functions 

/»oo 

u(x)=x, v(x)= (p t (0) ~pt(x))dt, w(x) = E x \X(t d )\ 

Jo 

are regular harmonic in D\ for u this is just the martingale property of X t , for v (the 
compensated potential kernel of X t ) this is proved, for example, in (32], while for w it 
follows directly from the definition. Furthermore, w(0) = u(0) = w(0) = 0 and v(x) = 
v(—x), w(x) = w(—x). It is known that 

v(x) ~ c 3 min(|x|, It^ 1 ) 

for x G R, with c 3 = 03 ( 01 , C 2 , ct) (see, for example, Lemma 2.14 in |19] ). I 11 particular, 
v(x) ~ c 3 |t| for x G D. Finally, by the boundary Harnack inequality given in Theorem 12. II 
(see Examples 5.5 and 5.13 in [12] for a detailed discussion), we have 

w(x) ~ C 4 v{x) ~ C 3 C 4 |a; | 

for x G (— -, |), with 04 = 04 ( 01 , 02 , a). Let us define 

f(x) = w{x) +u(x) = 2E x (|X(t d )|1[i i00 )(X(tz)))), 
g(x) = w(x) - u(x) = 2E a .(|X(r D )|l ( _ 00 _i](X(r D ))). 

Then / and g are non-negative, regular harmonic in D and equal to zero in (—1,1) \ D — 
{0}, so that they satisfy the assumptions of Theorem 13.11 On the other hand, 

f(x) f(~x) w(x) + x w(x) — x 4 xiv(x) 

g(x) g(~x) w(x) — x w(x) + x (w(x)) 2 — x 2 

for x G D. Since t/{t 2 — x 2 ) is decreasing in t G (x, 00 ), and w(x) < c :i c 4 x for x G (0, |), 

we obtain 

f(x) _ f(-x) > 4c 3 c 4 

g{x) g(-x) ~ (c 3 c 4 ) 2 - 1 

for x G (0, |). In particular, the limit of f{x)/g{x) as x —» 0 does not exist. 

4. Proofs of main results 

In this section we prove Theorem 13.11 We will always assume that To, R and D are 
fixed, where io G L, 0 < 2 R < R 0 and D C B{x 0 , R) is an open set. It is also understood 
that To G dD, although, at least formally, the argument extends also to t E D and 
t ^ D. Recall that the notation f ~ eg stands for c~ 1 g < f < eg with c > 0. 

We denote B r = B(x 0 ,r), B rs = B s \ B r , D r = D D B r and D rs = D s \D r when 

0 < r < s < R. We furthermore define D roo = D r r U (j£ \ Br). For a non-negative 

function / we let 

M r ,oo{f)=[ f{y)v(x 0 ,y)m(dy), M r>a (f) = [ f(y)u(x 0 ,y)m(dy). 

JX\B r J Br jS 

Finally, we let sd(x) = E x r D . 

To simplify the notation, we drop D from the notation in subscripts whenever possible, 
and we write r r = r Dr , t T)S = r Dra , s r (T) = s Dr (x), l r , s (T) = I^^t) etc. 
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Our argument is based on the boundary Harnack inequality of pj2], stated in Theo¬ 
rem o Under the assumptions of Theorem 13.11 the constant Cbhi^O; r, 2r, 3r, 4r) can 
be chosen so that it does not depend on r, as long as 0 < 4r < R, and it will be de¬ 
noted simply by Cbhi (recall that Xo and R are fixed). In a similar way, we denote 
0 L evy = C-uyy(x 0 ,r,2r) (with 0 < 2r < R 0 ) and C L evy-int = C h ^ y . int (x 0 ,r,2r) (with 
0 < 2r < Rq), chosen independently of r. With one exception, we will only use constants 
Cbhi, CLevy and Ci^vy-mt with these parameters. 

We prove Theorem 13.11 by considering separately two types of boundary points, which 
are called accessible and inaccessible in HU. First, however, we introduce some further 
notation and prove preliminary estimates. 


4.1. Decomposition of harmonic functions. From now on / and g are functions 
satisfying the assumptions of Theorem 13.11 and we assume that neither / nor g is equal 
to zero almost everywhere. Note that this implies that / and g are strictly positive in D. 
Whenever 0 < r < s < R, we decompose / into the sum of two functions, f TyS and / r>s , 
which correspond to the process X t exiting D r near its boundary (into D rtS ) and away of 
its boundary (into D s oo ): 


fr,s{x) = E x ((fl r , s )(X(r r ))), f r , s (x) = E x ((/l s ,oo)(X(r r ))). 

Not unexpectedly, a similar notation is used for the function g. Clearly, / = f T}S + f T)S , 
and both f rtS and / r;S are non-negative regular harmonic harmonic in D r which are equal 
to zero in B r \ D r . Therefore, we can apply Theorem 12.11 to f\ T)S and whenever 
0 < 4 r < s < R. 

Note that by Theorem 12.11 (with r — j), we have 

f(x) ~ CbhI-^3R/4,oo (/)E xT 2R/A 


for x G D r / 4 . Therefore, 


Mr }S (f ) ~ CbhI-^3H/4,oo(/)-^t,s( s R/2) (4.1) 

whenever 0 < r < s < j. The next result states, in particular, that there is little 
difference whether we write s R / 2 or s R in the above estimate. 

Lemma 4.1. If 0 < 8r < R, then 

E x T4 r U Ej.Tgy. ^ (1 T CBHlCLevyCj j g V y_j n t)E 3; T4 r . 


for x G D r . 

Proof. The hrst inequality is clear. For the other one, we use strong Markov property 
and Theorem 12.11 

Ea;Tg r E x T4 r E^Sgj. (N (T4 r )) 


< CbhiEzT^t- 

/ E y s 8 r(X(T 4 r))is(xo,y)m(dy) 


Jx\B 3r 

< CBHlE x T 4r 

/ E y T 8r v(x 0 ,y)m(dy). 


JX\B 2r 
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Furthermore, by Proposition 2.1 in [PS] (combined with the last displayed formula in the 
proof of this result), 


' X\B 2r 


E y T 8r v(x 0 ,y)m(dy) < ( sup E x t r . 


xGX 


>x\b 2 r 


u (xq, y)m(dy) 


— C^Levy 


fx\B 2r iy ( x 0 ’ y) m ( d y) 


Ix\B 16r u ( x o^y) m (dy)' 


It remains to use m- 


□ 


For convenience, we denote 

C T = 1 + C , BHlC'LevyC , Levy-inO 

so that S 4 r (x) ~ C T s 8r (x) if 0 < 8 r < R and x G D r . 

Our next result compares f 8r ,s with f 8r ,s- For f 8r ,s, we will use Theorem 12.11 which 
states that in D 2r we have f 8 r, s ~ CBmM erj00 (f 8rtS )E x T 4 r . The same estimate can be 
written down for f 8l . jS . However, M 6rj00 (f 8rjS ) involves an integral of f 8rs over -D 6r8r , 
which is often problematic. A much better estimate for f 8r ^ s can be easily obtained from 
the following corollary of Dynkin’s formula for X t . 


Lemma 4.2 (formula (2.12) in [12]). Let D C X be open and bounded, and let f be a 
non-negative function equal to zero in D. Then 

E J(X(r D )) = E :c r [ u{X t ,y)f{y)m{dy)dt (4.2) 

JO JX\D 

for x G D. 

Using the definition of f 8r)S and (12.11) to substitute u(x 0 ,y) for u(X t ,y) in (14.21) . we 
have 

fsr,s( x ) ~ CuevyOU), Sr, s)M S}00 (f )E x T 8r . (4.3) 

Note that not only we have M soo (/) instead of M 6T . ]00 (/ 8r . jS ), but also the constant 
FAevy (x 0 , 8 r, s) tends to 1 as r —^ 0 + . 

Lemma 4.3. IfO < 8r < s < j, then 

f8r,s(%) ^ Mor,s{.SR/f) 

f8r,s{' X ) ~ BHI 1 + M S ,R/a{s R/2 ) 
for x G D 2r . If 0 < 16r < s < then 

f8r,s{x ) ^ ^_3 r y_i r ,_ 3 XI 8rs (SiJ/ 2 ) 

/*,.(*) - BHI L6vy r 1 + M s , r/ 4 (s r/2 ) 

for x G D r . 

Proof. By Theorem 12.11 

/8r,s(-J') A C , BHI-^p6r,oo(/8r,s)E x T4 r , 

/'8r,s(j') ^ ^BHI'^^6 r ', 00 ( f8r,s)’^x^~4r m 

Furthermore, 

^Ior,oo\f8r,s) ^do r , S \f&r,s) A Mq r,s\f)i 

Mor,oo(f8r,s) A Ms,oo{f8r,s) = Ms,oo{f ) A ^3R/4,oo(/) + Af s ,R/4,(f )■ 
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Finally, by (14.11) . 

Mq r ,s(f) < CBmM 3 R/4 !00 (f)M er ^ s (sR/2), 

M s ,R/4:(f) > C B ft 1 M 3R /4 t00 (f)M St R/4(sR/2). 

We conclude that 

fsr,s{ x ) ^ ^-,4 -^6r,s( s R/2) 

fsr,s( X ) ~ BHI ! + Ms,R/4(s R / 2 ) ’ 

which is the desired upper bound. The lower bound is proved in a somewhat more 
complicated way. By Theorem 12.11 and estimate (14. 3 p , 

f8r,s( x ) ^ C'bhi M(ir ,oo ( fsr ., s ) 1~4r 5 

/ 8 r,s(^) < ^LevyM* ,00 (/)E x T 8r 

(we can write CWy = CLevy( 24 n 8 r, 16r) in the second inequality because s > 16r). By 
Lemma 14.11 E x T 8r < C r E x T 4 r . Furthermore, by Theorem 12.11 (as in (14.11) . but with R 
replaced by R/ 3) and again Lemma 14.11 

^6r,oo(f8r,s) = -^6r,s(f8r,s) ^ ^8r,s (f8r,s) = M 3r ,s(f) 

— C , BHI J ^fi/4,oo(/)^8r,s(Sfi/6) 

— C' B n l C T 2 MR/ 4::00 (f)M Sr ^(s2R/ 3 ). 

On the other hand, by (14.ip . 

M St00 (f) = M a , R/A (f ) + M r/4 ,oc(/) 

< C'bhI-^3R/4,oo(/)-^s,R/4(sr/2) + Mr/4, 00 (f) 

< M fl/4l0 o(/)(l + C Y BfflTf Si _R/4(Sij/2)). 

We conclude that 

f8r,s(%) ^ ri—3 /m—1 /^-3 ^8r,s( £ R/2 ) 

- BHI T 1 + m», k/4 ( Sh/2 ) ’ 

as desired. □ 


4.2. Inaccessible boundary points. Throughout this part we assume that xo is mac- 
cessible , that is, 

M 0 ,ooM = / E y T fl z/(x 0 ,2/)m(d2/) < oo. 

Jd r 

In this case /g riS and gsr,s turn out to be negligible compared to /sr-,s and gsr,s for suffi¬ 
ciently small r and s. 

Clearly, M 0:OO (sr/ 2 ) < M 0 )OO (sr) < oo. We remark that by (14.IF 

Mo,00(f) = M 0 , R/ 4 (f) + M r/4)00 (/) 

< OBHlTf3i?/4 )0 o(/)foro,R/4(s_R/2) + Mr/ 4 j0Q (/) < OO, 

and M(i^(g) < 00 by the same argument, and hence one can pass to the limit separately 
in the numerator and the denominator of (13.21) . 

Let 0 < e < 1. By the upper bound in Lemma [4.31 there is s = s(£) < eR such that if 
0 < 8 r < s, then 


fsr,s( X ) < £ fsr,s ( x ), 


g8r,s( X ) < £ 98r,s( x ) 


( 4 . 4 ) 
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for x G D 2r - Furthermore, estimate (14.31) and assumption lim Ci^vy^o , r , R) = 1 imply 
that there is r = r{e) < s/8 such that r_>0+ 

fsr,s( x ) ~ (! + e)E x T 8r M St00 (f), gsr,s ( x ) ~ (1 + e)E x T Sr M St00 (g ) (4.5) 


for x G D Sr . It follows that 

f( x ) 


< (1 + £ )/«-,w <(i +£) 


gsr,s ( x ) ~ (! + £)E x T 8r M Si00 (5f) 
3^,00 (/) 


5 f(x) gSr,s (j') M s ,oo(g ) 

for x G F ) 2 r- The lower bound is proved in a similar manner, and we obtain 


(1 + e) 


_ 3 M Si00 (/) /Or) ^ n , ^ 3 M ai0 o(/) 


< (l + e) J 


M Sj00 (g) g(x) M Si00 (g) 

for x G D 2x • Since £ was arbitrary and s converges to 0 as £ —» 0 + , we have 


(4.6) 


lim M = lim 

x ^ x -0 g[x) 


xGD s_>0+ M Sj00 (g) 

and Theorem 13.11 for inaccessible boundary points is proved. 


4.3. Accessible boundary points. In the second part of the proof we assume that Xq 
is accessible, that is, 

M 0:OO (s R ) = / E y r R u(x 0 ,y)m(dy) = oo. 

Jd r 

In this case f$ r , s and g% r .s dominate f% r)S and g$ r , s for all sufficiently small r. 

We remark that by (14.ip and Lemma 14.11 

M 0 ,oo(f) > Af 0 ,R/ 4 (/) > C B ^ 1 M 3R /^ oo (f)M 0yR / 4 (s R / 2 ) = oo, 

and Mo )OC (g) = oo by the same argument. In other words, the numerator and the 
denominator of the right-hand side of (j3.2[) diverge to infinity as r —» 0 + . In particular, 
if the limit of f(x)/g(x) in (13.21) exists, then it is automatically equal to the right-hand 
side. 

Onr argument is based on the following standard oscillation reduction lemma. 


Lemma 4.4. //'() < 8r < s < R 0 , then 

f&r,s(y) 


sup — inf 

y&D 2r y &D ^ 


< 




g»rM ^BHI T ^ 


sup — inf 

y£D s y^ D - 


Proof. For simplicity, we denote 

, f(y) 

A = sup —— , 
y&D s g{y) 

f(y ) 


B = sup 


f(y ) 
g{y)' 

f8r,s{y ) 


a = inf 


v* D * g(y) ’ 


b = inf 


^ 4^ / \ 1 

y&D 2r g%r,s\y) 

f8r,s{y ) 


y&D 2r gsr,s{y ) 


Since 


we clearly have 


(ig^-8r,s — f^-8r,s — Ag\g r ,si 


f8r,s — Ag^ rs . 


( 4 . 7 ) 
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In particular, a < b < B < A, and Theorem 12.11 applies to everywhere non-negative 
functions f 8rs — ag 8r ^ s , Ag 8r . jS — / 8r . )S and g 8rjS (note that / — ag and Ag — / typically fail 
to be non-negative everywhere). By 02.41) . 


f8r,s(y ) a 98r,s(y) 

sup - 7 T- 

y€D 2r 98r,s\y) 

A98r,s(y ) f8r,s(y ) 

sup - 7 T- 

y£D 2r 98r,s(y) 


< r 4 

— '~'BHI 

< c 4 

— °BHI 


inf 

y&D 2r 


inf 

y€D 2r 


f8r,s(y ) a 98r,s(y) 

98r,s{y ) 

^98r,s{y) f8r,s(y ) 
98r,s(y ) 


This translates to B — a < C^ m {b — a ) an d A — b < C^ m (A — B), and adding the sides 
of these inequalities leads to the desired inequality 

(^Bffl + 1 )(-B — 5) < (CgHi — 1)(^4 — a). □ 


For continuous processes (in sufficiently regular domains), the above lemma easily yields 
the assertion of Theorem 13 .1 1 For processes with jumps one needs to incorporate the non¬ 
local parts f 8riS and g 8rtS using Lemma [4.31 As it was remarked in the introduction, this 
modification was developed in |j7], and extended in HU- 

Let 0 < e < 1 and 0 < s < ^|. By the lower bound in Lemma 14.31 there is r = r(e, s ) < 
| such that 

o 

f8r, S ( X ) < ^/Sr, s (^), 98r,s( X ) < £ 98r,s ( X ) ( 4 -8) 


for x G D r . It follows that 


sup — inf 

x£D r xeD r 


f( x ) . n . X fsr,a( X ) 1 • f fsrA X ) 

—— < (1 + e) sup —- -- mf —V — 

9\x) xeD r 98r,s\ x ) 1 + £ g 8 rA X ) 


By Lemma [4.41 and the inequality 1 — (1 + e) 1 < e, 
f( x ) ^ ^bhi — 1 • r \ f ( x ) 


sup 

XdD r 


inf 

xGD, 


< 


g( x ) 0 BHI + 1 


sup 

x&D s 


inf 

xGD. 


g( x ) 


+ e ( sup + inf 


XdD r 


x£D 7 


f 8 r,s ( x ) 


98r 


[X 


(4.9) 


Denote by Q the upper limit of the expression in the left-hand side as r —> 0 + . Using (14.71) 
and taking the upper limit of both sides as s —> 0 + leads to 


Q < ((° H1 — Q + 2 e sup 7M 


^BHI + 1 


that is, 


Q <e(l + C BHI ) sup 


xeDfl/4 9\ x ) 

f( x ) 


x&d r/ 4 g( x ) 

Since £ is arbitrary, we conclude that Q = 0, and the proof of Theorem 13.11 is complete. 


4.4. Extensions. We first prove the statement contained in Remark 13.31 Denote g(x) = 
E x tr. Then g is not a regular harmonic function in Dr, but for every open U C Dr, 

g(x) = E x Tu + E x g(X(Tu)). 

We interpret E x tr as if it originated from a jump to a distant point (a point at infinity), 
and we define 

Mr,oo(g) = 1 + / g(y)u(x 0 ,y)m(dy), g r , s ( x ) = E x r r + E x ((gl Sj00 ){X(r r )))-, 

Jx\B r 

the definitions of M r s (g) and g r . s { x ) for finite s remain unaltered. One can then follow 
carefully the proof of Theorem 13.11 and see that no changes are required. This shows 
validity of Remark 13.31 
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In the remaining part of this section we argue that an extension stated in Remark 13.41 
is true: the limit in Theorem 13.11 converges uniformly in / and g , and also in D, in the 
sense of (13.31) . 

We claim that if 0 < q < Ro and p > 0, then there is p , which depends only on q, p 
and the characteristics of the process X t , such that 0 < p < q and 


s u PxeD p (f(x)/g(x)) 1 f swp xeDq (f(x)/g(x)) \ 

inf xeD p (f(x)/g(x)) ~ C^ m +1 \M xeDq (f(x)/g(x)) ) 


for all open sets D and all functions / and g as in Theorem 13.11 (this estimate is very 
similar to (14.91) 1. By considering the supremum of both sides of (14.101) over all /, g and 
D , and then taking the upper limit as q —> 0 + , we obtain the desired result: 


lim sup sup 

r—>0+ D,f,g 


( sup xeDr (f( x )/9(z)) A 

V inf xeD r (f(x)/g(x)) ) 


< 


?/(l + Cbhi) 
2 


for arbitrary p > 0. Therefore, it remains to prove (14.101) . 

Let 0 < q < yj-Ro and p > 0. We consider two additional parameters 6, N > 0; 
the actual values of 6 (small real) and N (large integer) are to be specified at the end 
of the argument. By the assumption lim r _> 0 + CLevy(+ 0 j r, R) = 1 one can construct a 
decreasing sequence of radii ao,a ±,... ,ajv so that ao is the input radius q, |will be 
the output radius p, and we have 16a„+i < a n and (Wevy^'o, 8a n+] , a n ) < 1 + 5 for all 
n — 0,1,..., N — 1. 

Following HU, we consider two scenarios. Suppose first that for some n we have 

M an+1 ,a n (sR/2) < <^(1 + M antR / 4 (s R / 2 ))■ (4-H) 

Then the argument is fairly simple: as in Section 14.21 by Lemma l4~3l we have the inequal¬ 
ity (14.41) with r = a n+1 , s = a n and e = C^ m 5. Since 8 a n+ i, a n ) <1 + 5, the 

estimate (14.51) holds with r = a n+1 , s = a n and e = S. This implies (14.61) (with s = a n , 
x £ -D 2a „ + i and £ = Cb HI 5), and in particular the left-hand side of (14.101) does not exceed 
(1 +C'b HI 5 ) 6 —1. Estimate (14.101) follows with p = a n+ 1 , provided that (1 +Cb HI 5 ) 6 —1 < p. 
We choose S small enough, so that this inequality is satisfied. 

In the other scenario, for each n the converse of (14.111) holds. Summing up these 
inequalities for n = 0,1,..., N — 1 we obtain 


Ma N ,a 0 (s R / 2 ) > N5( 1 + M a0)R / 4 (s R / 2 )), 

and we argue as in Section 14.31 Again by Lemma 14.31 we have (14.81) with r = |cqv, 
s = ao and e = C , B HI C'Lev y C , r(^^) _1 - Inequality (14.91) follows. Dividing both sides of it 
by inf x eD r (f (x)/g(x)) and using monotonicity of this expression in r, we obtain (14.101) 
for p = 1 oat, provided that £(Cbhi + 1) < p- Since 5 is now fixed, we may choose N 
large enough, so that this condition is satisfied. This completes the proof of the extension 
described in Remark 13.41 


4.5. Martin representation. In this section we prove Theorem 13.61 We assume that 
the assumptions of Theorem 13. II are satisfied in a uniform way for all xq G D. 

We note one important property of the Green function: if U is an open subset of D, 
then 

G D (x,y) = E x G D (X(Tu),y) + Gu(x,y ) (4.12) 

(where, as usual, we assume that Gjj(x,y ) = 0 whenever x ^ U or y ^ U). In particular, 
Gr>{x,y) is a regular harmonic function in D \ B(y,r) for every r > 0. By a duality 
argument, Gd{x,u) is a regular co-harmonic function in D \ B(x,r ) for every r > 0. 
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Proof of Theorem l M a )[ The assumptions are completely symmetric under duality, and 
hence we may apply Theorem 13.11 to both harmonic and co-harmonic functions. In par¬ 
ticular, as already remarked before the statement of Theorem 13.61 the Martin kernel, 
defined as the boundary limit of co-harmonic functions 


M d (x,z ) 


UmkdM), 

*?ZG D (x,y) 


exists for all boundary points z G D (here and below x E D is a fixed reference point). 
In other words, the Martin boundary coincides with the Euclidean boundary. □ 


The representation given in part (d) essentially follows now from the general theory 
of Martin boundary, together with some ideas developed in m- For simplicity, in the 
remaining part of the proof we simply write that a function is harmonic when we refer 
to harmonicity in D. 


Proof of Theorem! 3. C\b) . Following the proof of Theorem 2 in [[IT], we find that 


M d (x , Xq ) is a harmonic function with respect to x if and only if x 0 is accessible. Indeed, 
for an inaccessible boundary point Xq we have, by (13.21) in Theorem 13.11 


M d (x,x 0 ) = C u(y,x 0 )G D (x,y)m(dy) 

J D 

for C = (f D u(y, x 0 )Gd(x, y)m{dy))~ l > 0, and so the Martin kernel is not harmonic (to 
see this, simply use (I4.12D and Fubini). On the other hand, if Xq is accessible and R > 0, 
then 

G d(X (t d \b( xo ,r)) > v) 


E x M d {X{t DX b {xo>r) ),xq) = E, 1 i m 

yeD 


G D (x,y) 


(4.13) 


Recall that Gd(x, y ) is a regular harmonic function of x 6 D\B{x o, R) when y G B{x o, R). 
By Fatou’s lemma, 


^ x M d (X(t d ^ xo R) ),xq) < M d (x,x o), (4.14) 

and we claim that in fact equality holds, that is, we can exchange the limit with the 
expectation in (14.131) . By Vitali’s convergence theorem, it suffices to prove that that the 
ratio in the right-hand side of (14.131) is a uniformly integrable family of random variables 
for y G D D B(x o, r) for some r > 0. The argument is exactly the same as in the proof of 
formula (77) in |11| ; for the convenience of the reader, we repeat it below. 

Assume that 0 < 8 r < R and that x, x ^ D D B(x o, R). We will first prove that 


Gd(z, V ) 
sup n (r v 

y£DnB(xo,r) '^ r D\X, y) 
z£D\B(x o,4r) 


< oo. 


(4.15) 


By the boundary Harnack inequality (Theorem 12 .11) applied to Gp(z, •) and Gd(x, •), it 
suffices to consider a fixed y G D D B(x 0 ,r), that is, to show that Gn{-,y) is bounded 
in D \ R(xo,4r). This is relatively simple, but somewhat technical. Denote Di = D D 
B(x 0 ,r), D 2 = D fl B(x 0 ,r), D 4 = D D B(x 0 ,Ar) and D r = D\ B(x 0 ,4r). By Dynkin’s 
formula fl4.2[) . 


E z (G j d(X(td/), y)lD 2 (X(rD'))) 


< ( sup v{v, w) 

D' 

WGD2 


ID' 


Go'(z,v)m(dv) / Gr)(w,y)m(dw). 

Jd 2 
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The supremum is finite by Assumption O and boundedness of D , and the integrals in 
the right-hand side are bounded by sup. ueD E u To and sup, ieD E u td, respectively. Further¬ 
more, 

^z{GD{X{T D i) 1 y)l Di \ D2 {X[r D t))) < sup Gd(v,w), 

v£D4\D2 

w£Di 

and the right-hand side is finite by Assumption [0 By adding the sides of these two 
bounds and using harmonicity of the Green function, we complete the proof of (14.151) . 

On the other hand, if we denote D" — D \ B(xq,8t ) and D'" — D \ Bfxo,R), then, 
again by Lemma 14.21 

E x (Gd (X ( td ), y) 1_d 4 (X (t£>/// ))) 


A C\.£ V y 


p(t 0 , v)Gd'"(x, v)m(dv) 


ID"’ 


Gr,{w , y)m(dw) 


iDi 


and, in a similar way, 

Go(x,y) > (Gd(A(t£,»), y)l D4 (X(Tr>"))) 


> C£ y 


z/(xo, v)Gd"(x, v)m(dv) 


ID" 


G D (w,y)m(dw] 


id 4 


It follows that 

E 


G D(X(T D \ B{xoJl) ),y) 

- G ^~ y j - i-DnB(xoAr){^{ T D\B(x 0 ,R))) 


<cl u(x 0 ,v)G D \B( X0 ,8r)(x,v)m(dv) 

\J D\B(x o,8r) 

where C does not depend on (sufficiently small) r > 0 and y <E D D B(x o, r). Recall that 
Gd\b( x 0 ,8r)(x,v) increases to Gjj(x,v) (because the corresponding exit times td\b(x 0 ,8t) 
increase to To). By monotone convergence, the right-hand side converges to zero as 
r —» 0 + . Together with (14.151) . this completes the proof of uniform integrability of the 
right-hand side of (14.131) . 

Part (b) follows, and in addition we see that for accessible boundary points z, the 
Martin kernel Mo{x, z ) is a regular harmonic function in D\B(z, r ) for every r > 0. □ 


In order to apply the general theory of Martin boundary, we need to prove 
that the Green operator, which maps a measurable function f(x) to Gof(x) = 
J D Gd(x, y) f (y)m(dy) , takes bounded functions into continuous ones. Let / be a bounded 
function on D, xq e D and e > 0. Clearly, \Gof(x)\ < ||/||E x td for x G D, so that Gof 
is bounded. Let r > 0 be small enough, so that E x ' T B{xo,r) < £ for x G B(x 0 , r ). By (14. 12[) . 

GdJ(x) = Fi x Gof {X(t B ( xqir ))) + GB(xo,r)f ( x ) ■ 

The first term is continuous in B(xq,t) by Theorem 13.II (see Remark 13.51) . The other one 
is bounded by e||/||, an arbitrarily small number. Therefore, Gof is continuous at Xq. 

The general theory of Martin boundary tells us now that if / satisfies the assumptions 
of Theorem 13.61 and / is equal to zero in the complement of I), then 

f(x)= f M D {x,z)ii(dz) (4.16) 

JdmD 

for some measure /i on the set of accessible boundary points d m D, see Theorem 14.8 
in |n>]. Furthermore, if we show that for every z G d m D, Mo{x, z) is a minimal harmonic 
function with respect to x, then the measure /i in the above representation is unique. 
Minimality of Mo{x, z) is proved as in the final part of the proof of Lemma 14 in m- 
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Proof of Theorem \3.tf(c)\ Suppose that / is harmonic, 0 < f{x) < M D (x, x 0 ) for all 
(in particular, f(x) = 0 for x G 3C\D) and that the measure y in representation (14.161) 
is zero on d m D fl H(xo,4r) for some r > 0. Our goal is to prove that / is identically 
zero. This will imply that if / is harmonic and 0 < f(x) < Md(x,x o) for all i£l, then 
the measure y in representation (14.161) is concentrated in {xo}, and thus Md(x,x o) is a 
minimal harmonic function. 

For every z G d m D\B(x o, 4r), Md(x, z ) is a regular harmonic function in DP\B(x o, Sr). 
Hence, by Fubini, / also has this property. Furthermore, by the boundary Harnack 
inequality (Theorem 12. ip . / is bounded on D 0 B(x 0 , 2r). 

On the other hand, since f(x) < Md(x,x 0 ), one easily hnds that / is also a regular 
harmonic function in D \ B(xo,r). This is exactly the same argument as in Lemma 9 
in DU; for the convenience of the reader, we provide the details at the end of this section. 
In particular, since / is bounded in D fl B(x o, 2r), it is bounded on D. 

A sweeping argument, which is a simplified version of Lemma 10 in DU, proves then 
that / is a regular harmonic function in D\ Let cr n be the sequence of consecutive exit 
times from alternately D fl B(xo,4r) and D \ B(xo,r). That is, cr 0 = 0 and a n+ \ = 
cr n + T v ° $<T n , where V — D fl B(x o, 4r) when n is even and V = D \ B(x o, r) when n is 
odd (and is the shift operator). 

Clearly, cr n < To < oo. Since a n is increasing, by quasi-left continuity, X(cr n ) has a 
limit as n —> oo. Therefore, it is impossible that cr n < T£, for infinitely many n. It follows 
that with probability one, eventually a n — Td- 

Since f{x) = E x f(X(a n )) and / is bounded, by dominated convergence we have f(x) = 
E x f(X(jD)) = 0, as desired. □ 


We have thus proved the representation (14.161) for harmonic functions / which are zero 
in the complement of D. The general case is handled as in Lemma 13 in HP- 


Proof of Theorem \3jL |(d) , Let D n be an ascending sequence of open sets such that D n C D 
and U£U D n = D. Then, by Lemma [4721 


f(x) = E x f(X(r Dn )) > [ ( [ G Dn (x,y)u(y,z)m(dy)) f(z)m(dz). 

JX\D \JD n J 

The integrand in the right-hand side increases as n —> oo, and therefore by monotone 
convergence, 


f(x)> j^ (^J G D {x,y)v{y,z)m{dy)^ f{z)m{dz). (4.17) 

Let g(x) be equal to the right-hand side of (14.17)) for x G D, and to f{x) for iGl\h. 
From Lemma 14.21 and the property (14.121) of the Green function it follows easily that g is 
harmonic: if U is open and 1 /C£), then 


Exg(X(Tu)) = 'E, x (fl x \ D )(X(T U )) + E X [ (! G D (X(Tu),y)u(y, z)m(dy)\ f(z)m{dz) 

J X\D \Jd J 

= / ( [ (Gu(x,y) + E x G D {X(r u ),y))u(y,z)m(dy)\f(z)m(dz) = g(x). 

J X\D \J D J 

Therefore, f —g is a non-negative harmonic function which is equal to zero in X.\D, and 
so it has a unique representation (14.161) . 

Finally, the outer integral in (14.171) is finite, and so points at which the inner integral is 
infinite cannot contribute to the integral. It follows that we can change the outer integral 
to an integral over X\(D U d m D). The proof of (13.41) is complete. □ 
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Proof of Theorem 1 .7. c) [ By the boundary Harnack inequality, if the right-hand side 
of (13.41) is finite at some x E D, it is finite everywhere in D. Indeed, let / be given 
by (13.41) . If f(x) = oo for some x E D, by Theorem 12.11 f is infinite at every point of a 
ball B(x,r ) contained in I). If y E I), then again using Theorem 12. 1 1 (for a ball centred 
at y), f is infinite at y. 

Finally, harmonicity of the right-hand side of (13.41) . whenever it is finite, follows from 
Lemma 14.2[ property (14.121) of the Green function, and harmonicity of the Martin kernel. 

□ 

At the end of this section, we present the proof of Lemma 9 in HU, adapted to our 
setting. This result was used in the proof of Theorem I3.£|lc)| 

Lemma 4.5 (Lemma 9 in [TT|). Let U and D be open subsets of X such that U C D. 
If 0 < f(x) < g(x) for all x E 3L, f and g are harm.onic in D, g is a regular harmonic 
function in U and g(x) = 0 for x E £ \ D, then f is a regular harmonic function in U. 

Proof. Let D n be an ascending sequence of open sets such that D n C D and Ll Dn = D, 
and let U n = U D D n . Then rjj n increases to tjj, and, by quasi-left continuity, X{ju n ) 
converges to X(tu) with probability one. It follows that if X(tu) E D\U, then eventually 
Tu n = tjj for n large enough up to an event of probability zero. Hence, 

lim E x (gl D \ u )(X(r Uri )) = E x (gl D \ u )(X(r u )) = g(x). 

n—>oo 

Therefore, 

E x (fl UXUn )(X(r Un )) < E x (glu\u n ) (X(r Un )) = g(x) - E x (gl D \ u )X(r Uri )) 
converges to zero as n —> cxd. It follows that 

f(x) = lim E x f(X(T Un )) = lim E x (fl D \u)(X(T Un )) 

oo n—> oo 

= E x .(/1 dw )(X(t c/ )) = E J(X(tu)), 

as desired. □ 
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